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ABSTRACT 

We investigate the hypothesis that the cores of cUiptical galaxies and bulges are created 
from the binding energy liberated by the coalescence of supermassive binary black 
holes during galaxy mergers. Assuming that the central density profiles of galaxies 



were initially steep power laws, p 



we define the "mass deficit" as the mass in 



o 



stars that had to be removed from the nucleus in order to produce the observed core. 
We use nonparametric deprojection to compute the mass deficit in a sample of 35 
early-type galaxies with high-resolution imaging data. We find that the mass deficit 
correlates well with the mass of the nuclear black hole, consistent with the predictions 
of merger models. We argue that cores in halos of non-interacting dark matter particles 
should be comparable in size to those observed in the stars. 

Key words: black holes: binary black holes — galaxies: elliptical and lenticular, cD — 
galaxies: interactions — galaxies: nuclei 



INTRODUCTION 



ing energy to the surrounding stars ( [ jjegelman, Blandford 



Ferrarese et al. (1994) and |Lauer et al. {1991 ) divide ellipti 



Rees 198C 



Ebisuzaki, Makino & Okumura 1991). High- 



cal galaxies into two classes based on their nuclear proper- 
ties, which Lauer et al. call "core" and "power-law" galaxies. 
Core galaxies exhibit a definite break in the surface bright- 
ness profile at some radius Ri,; inward of this break, the 
logarithmic slope gently decreases in a manner that mim- 
ics a constant-density core. Power-law galaxies show essen- 
tially a single power-law profile throughout their inner re- 
gions, T.{R) ~ R-^, r ^ -0.8 ±0.2. The brightest galaxies, 
Mv ^ —21, are exclusively core galaxies while faint galaxies, 
Mv ^ — 16 always exhibit power laws; galaxies of intermedi- 



resolution A^-body simulations verify that this process can 
convert a steep power-law cusp, p ~ r~^, into a shallow 
power-law cusp, p ~ r ~^, within the radius of gravit ational 
influence of the BHs (iMilosavljevic & Merritt 200l|). Suc- 



ate lun inosity can exhibit either type of profile ( [CJcbhardt 
et al. lyyb T- wnue tne two categories were mitially seen 



cessive mergers would presumably lower the density of the 
core still more. In this model, power-law galaxies are those 
which have not experienced a major merger since the era of 
peak BH growt h, or which have re-generated t heir cusps via 
star formation (Milosavljevic & Merritt 2001). 

Preli minary tests of the cu sp-di sruption model were pre- 
sented by Faber et al. (1997) and Milosavljevic & Merritt 
(2001). The former authors plotted core properties (break 



as distinct, nonparametric deprojection revealed that even 
the "core" galaxies exhibit power laws in their central spac e 
densities, p ~ r"^, with 7 1 (Mer ritt fc Fridman 1995 ). 



Power-law galaxies have 1 < 7 < 2.5 ( [Cebhardt et al. 1996| ). 
Furthermore the distribution of de-projected slopes is essen- 
tially continuous as a functi on of galaxy luminosity in the 



radius, core luminosity) versus global properties in a sample 
of 19 early-type galaxies and noted a rough proportionality. 
In the paradigm investigated here, core properties should 
correlate more fundamentally with BH mass, since the mass 
of stars ejected by a de caying binary BH is expected to be 



of order the BHs' mass. Milosavljevic & Merritt (2001) used 



larger samp les now available ( Ravindranath et al. 2001 



etaL2mj 



Rest 



the new empirical relation between gal axy velocity disper- 



sion a nd BH mass, the A/. — a relation ( Ferrarese & Merritt 
200c ; GJebhardt et al. 2000), to estimate BH masses in the 



Here we assume that the steep central density cusps of 
faint ellipticals and bulges, p ~ r~^, are characteristic of 
the earliest generation of galaxies, and ask: How do the low- 
density cores associated with bright galaxies form? An ap- 
pealing hypothesis links cores to nuclear black holes (BHs): 
in a galactic merger, the BHs will fall to the center of the 
merger remnant and form a bound pair, releasing their bind- 



Faber et al. sample. They found that a rough dynamical es- 
timate of the "mass deficit" - the mass that would need to 
be removed from an initially density cusp in order to 
produce the observed profile - correlated well with M,. 

In this paper we present the most careful test to date of 
the BH merger hypothesis for the formation of galaxy cores. 
Wc use nonparametric deprojection to compute the mass 
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deficit in a sample of galaxies with high-resolution imaging 
data from HST (§2). We find a strong correlation between 
this mass and the mass of the nuclear BH; typical ejected 
masses are ~ several M, (§3). We argue (§4) that this re- 
sult is consistent with the formation of cores via hierarchical 
mergers of galaxies containing pre-existing BHs. Cusps of 
non-interacting dark matter particles should behave in the 
same way as cusps of stars in response to heating by binary 
BHs, and we argue that the damage done to stellar cusps 
by this mechanism is a reasonable guide to the damage that 
would be done to dark matter cusps. 



2 DATA AND METHOD 

Our data set is drawn from a sample of 67 surface 
brightness profile s of early-type g alaxies observed with 
HST/WFPC2 by ^est et al. (200l| ), and three additional 
galaxies: NGC 4472 and 4473, observed with WFPCl by 
Ferrarese et al. (1994| ); and a WFPC2 F547M image of M87 
( Jordan et al. 2002 ). The Rest et al. sample was selected 



from tlje'seT'or'all early-type galaxies with radial velocities 



less than 3400 km s~ , absolute V-band magnitudes less 
than —18.5, and absolute galactic latitude exceeding 20 de- 
grees. From this sample we excluded 13 galaxies for which 
central velocity dispersions were not available in literature; 
as discussed below, velocity dispersions were needed to com- 
pute BH masses in most of the galaxies. For specifics of the 
image-data reduction we refer the reader to the sources cited 
above. Surface brightness profiles used in this study were 
major-axis profiles. We applied a crude correction for the ap- 
parent ellipticity of the galaxies by multiplying the volume- 
integrated quantities, defined below, by (1 — tb), where eb is 
the ellipticity of the isophote at the break radius. 

The intrinsic luminosity profiles i^(r) were obtained by 
deprojecting the PSF- deconvolved surface brightness pro- 



files S(_R) using the n on-parametric MPL technique (Mer- 



ritt & Tremblay 1994). We opted for one-step deprojection 



via maximization of the penalized likelihood functional 



- A 



poo 


(f log V 


'o 


(dlog r-)2 



ci log 



(1) 



The first term compares the observed surface brightness data 
Ei = T,{Ri) to the projections Pilv] of the intrinsic luminos- 
ity estimate given by the operator 



u{r)rdr 



(2) 



The second term in equation ^ is the penalty function which 
assigns zero penalty to any power-law !^(r); hence our esti- 
mate of 1/ is unbiassed if i/ is an unbroken power law and 
should be minimally biassed if i/ is approximately a power 
law. The relative strength of the penalty term is regulated 
through the parameter A which was chosen by eye and equal 
to 0.01 for all galaxies; integrated quantities like the mass 
deficit defined below are only weakly dependent on A. 

Distance s for 26 of the g alaxies were drawn from the 

SBF survey ( rnnry pt rd 9001 ) Fnr tVip rpmaining ."^4 galav- 
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Figure 1. Definition of the mass deficit, illustrated using the pro- 
file of NGC 5903. Left panel: density in Mqpc~^ as a function of 
radius in parsecs (thick line); hypothetical original density cusps 
for fiducial slopes 70 = 2, 1.75, 1.5 (thin lines). Dots indicate the 
corresponding break radii. Shaded region is the mass deficit for 
70 = 1.5. Right panel: the negative logarithmic derivative of the 
density profile, 7 = — dlog p/d log r. Break radii are defined as the 
radii where 7 = 70. 



the individual m ass-to-light ratios Tv quoted in Magor 



rian et al. (1998) or their best-fit relation log(Tv/T0,v) = 
-1.11 ± 0.33 + (0.18 ± 0.03) log(Lv-/iQ,v) for galaxies not 
included in that study. 

We define 7 = — dlog /9(r)/cilog(r) as the local, negative 
logarithmic slope of the deprojected density profile. Power- 
law galaxies are defined as those in which 7 > 2 at all radii; 
typically the profiles of such galaxies show no clear feature 
that can be identified as a "break radius" and are unlikely 
candidates for cusp destruction by binary BHs. In the re- 
maining 35, "core" galaxies, the slope varies from 7 > 2 
at large radii to 7 < 2 at small radii; the radius at which 
the slope crosses 7 = 2 in the positive sense {d'y/dr > 0) 
is called here the "break radius" rt (Figure ^). This defini- 
tion has little in common with the more standard definition 
based on fitting of the surface brightness profile to an ad 
hoc parametric function. In four galaxies the slope crosses 
7 = 2 in the positive sense at more than one radius and 
thus the definition of rb is ambiguous. In such cases, we se- 
lect the crossing toward larger radius from the largest dip of 
the slope below 7 = 2. 

We define the mass deficit as the difference in integrated 
mass between the deprojected density profile p{r) and a 7 = 
70 = 2 profile extrapolated inward from the break radius: 



Mdof = 47r(l - et,) 



- Pir) 



(3) 



Our choice of an r~^ density profile to characterize the 
"undisrupted" core is to a certain extent arbitrary; adia- 
batic growth of BHs can produce cusps with 1.5 7 2.5 
depending on initial conditions, and the faintest ellipticals 
with measured cusp sl opes exhibit a similar range of 7's (e.g., 
Gebhardt et al. (199(:)). To test the sensitivity of our results 
to the assumed initial profile, we repeated the analysis using 
fiducial slopes of 70 = 1.75 and 1.5. Values of 70 > 2 were 
found to exclude all but a few galaxies. 

BH masses for a few of the galaxies in our sampl e are 



available from spatial ly-resolved kinematical studies (Mer- 



les we adopted distances computed assuming a pure Hub- 



ritt & Ferrarese 2001). For all other galaxies we estimated 



ble expansion with Hp — 80 km s~^ Mpc~^ corrected for 
Virgo-centric infall (Rest et al. 2001). Luminosity densi- 
ties were converted to mass densities p(r) — Ti/(r) using 



M, via the Al, — a relation. 



M, ^ 1.4 



X 10* M0 



200 km s" 



(4) 
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( Ferrarese & Merritt 200C ) , where CTc is the central velocity 
dispersion corrected to an aperture of re/8 with re the effec- 



tive ra dius. We used the aperture corrections of j0rgensen 
et al. ( | L995) to co mpute ctc from published values of a in 
Davies eii ai. (1987] ), Tonry fc Davis (198l| ) and [Pi Nella~ei 



"(19^ 



3 RESULTS 



Results are given in Tables ^ and |^ and Figure ^; Table ^ 
gives mass deficits only for 70 = 2 while Figure ^ shows Mdof 
computed using all three values of the fiducial slope, 7o=(2, 
1.75, 1.5). In the first panel of Figure ^ we have also plotted 
dynamical" estimates of A/dcf for a sample of gala xies from 



Gebhard^_et_a/j_^199q) , using equation (41) from Milosavl- 
evic fc Merritt (2001|), Mdy„ = 2(2 - 7)/(3 - 'y)a^Rb/G 



Mdyn depends on the density profile only through Rb, the 
break radius of the surface bightness profile, and 7. When 
calculated for the deprojected galaxies in our sample, Mdyn 
was consistent within the scatter with Mdef- 

For 70 = 2, the mass deficits are clustered about a linear 
relation defined by (log(Mdef/M.)) = 0.92, 1.0, and 0.65, re- 
spectively, for Es and SOs; Es only; and SOs only, correspond- 
ing to Mdof ~ (8.4, 10, 4.5)Af, . Decreasing 70 decreases Mdcf 
(cf. Fig. lb) and Ad^d becomes negative/undefined in galax- 
ies when the minimum pointwise slope 7min approaches 70. 
We do not cite values of (log(i\fdcf/-Af,)) for these low values 
of 70 since the mean depends strongly on which galaxies are 
defined as having "cores." However Figures ^b,c shows that 
Mdef remains of order Al, or greater for many of the galax- 
ies even when 70 < 2. We speculate that the lower mean 
value of A/dof for the SOs may indicate a role for gaseous 
dissipation in the re-formation of cusps following mergers. 
We emphasize, however, that the mass deficit of lenticulars 
alone appears to be completely uncorrelated with the BH 
mass. 

Table |^ lists linear regression fits to the data that were 



carried out using the rou tine of Akritas & Bershady (1996) 
For the |Rest et al. (2001 ) ellipticals, as well as for the com- 



plete sample, the fitted power-law indices of the Mdcf — M. 
relation calculated with 70 = 2.0 and 1.75 are statistically 
consistent with unity. To test the sensitivity of the fitting 
parameters to the assumed power- law index a of the M. — a 
relation (equation W), we recomputed the fits with the Gleb- 



hardt et al. (200C ) value a — 3.75 ± 0.3 in addition to the 
Ferrarese fc Merritt (2000|) value a = 4.85 ±0.5 that is stan- 



dard throughout the current paper. The effect of changing 
to a shallower version of the M, — a relation is a steepening 
of the Mdef - M. relation, from d log A/dcf /d log M. « 0.91 
to d log Mdef /d log M. ~ 1-16. 

We note that in most galaxies the break radius defined 
via 7(ri,) = 70 = 2 is close to the radius at which the den- 
sity profile exhibits a visual break, i.e., where the curvature 
is greatest. The break radius, however, is not a good pre- 
dictor of the mass deficit; in particular, some galaxies with 
large break radii, rj, ^ r,, with r, = GM,/a^ the dynam- 
ical radius of the BH, have mass deficits below the mean. 
We speculate below that these large break radii may be pro- 
duced by mechanisms other than BH mergers. 

Although estimation of errors in Figure ^ is difficult, 
we believe that the scatter is at least partly intrinsic. Un- 
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Figure 2. Mass deficit vs. BH mass for three different values of 
70, the assumed logarithmic slope of the density cusp before en- 
ergy input from BHs. (a) 70 = 2; (b) 70 = 1.75; (c) 70 = 1.5. 
Filled circles are elliptical galaxies and empty circles are lenticu- 
lars. Squares show dyn amical estimates of Mam- for th e sample of 
galaxies considered in Milosavljevic & Merritt (2001). Solid line 
is the Mdot 



M, relation. Dotted line in panel (a) is a linear 
regression fit to all galaxies (see Table 3) . 
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certainties in log M, are due primarily to uncertainties in 
o" and are of order log 2. Uncertainties in log Mdef are also 
roughly log 2 based on variances between the redshift-based 
and SBF distances. A scatter greater than that due to mea- 
surement uncertainties would be reasonable given the diflter- 
ent merger histories of galaxies with given M, . Similarly, if 
the progenitors of the galaxies exhibited a range of different 
70, this could in itself explain the observed scatter in Figure 



4 DISCUSSION 



The mass ejected by a decaying BH binary is 



Moj « JMw In 



(5) 



where M12 = Mi + M2 is the binary mass, an is the semi- 
major axis when the binary first becomes hard, and the 
separation at which the rate of energy loss to grav itational 
radiation equals the rate of energy loss to the stars dQuinlan 



199(: ) . J is a dimens ionless mass-ejection rate; for e q ual-mass 
binaries, J « 0.5 (Milosavljevic & Merritt 2001). [Quinlan 



(1996) claims that J is nearly independent of M2/M1 even 



for extreme mass ratios, implying that a mass in stars of 
order Mi + M2 is ejected during every accretion event, even 
when M2 is tiny compared with M\ . This non-intuitive result 
is due to Quinlan's ejection criterion, which for M2 ^ Mi 
includes stars that would not have gained sufficient energy 
to escape from the binary. If instead we equate the change in 
energy of the binary with the energy carried away by stars 



that are ejected with v ^ Vun, we find Me 



M2. This 



argumelnt suggests a relation 



Moj f» M2 In 



ah 



(6) 



which is consistent with equation (5) when Mi « M2. 

Using equation (6), and adopting Merritt 's (2000) semi- 
analytic model for decay of a binary in a power-law cusp 
we find 

/Mi\«-2 - 



a 



gr 
ah 
and 



:A\lnA\' 



■ 4.6M2 



A« 7.5 



VM2y 



Ml < M2. 



(7) 



(8) 



Thus Mej/Ah varies only negligibly with M\/M2- Hence- 
forth we adopt Me-, « 5M2. 

If a BH grows by sequential accretion of smaller BHs, 
this result implies a mass deficit of order five times the final 
BH mass. However if the BH grows via a merger hierarchy 
of comparably-massive BHs, we expect Mde/ to be larger. 
The idea here is that the damage done to cusps is cumu- 
lative: a merger of two galaxies whose cusps had previously 
been destroyed by binary BHs, will produce a shallower pro- 
file than a merger between two galaxies with initially steep 
cusps, even if the final BH mass is the same. Galaxies with 
masses M ^ 10^^ A/0, including most of the galaxies plot- 
ted on Figure ^, are believed to have undergone at least one 
major merger since a redshift of 1 (e.g. Kauffmann, Chariot 
& Balogh 2001). Thus we predict Mdef ^ 5M., consistent 
with Figure [i] if 70 ^ 1.5. 



Our interpretation of the mass deficit depends critically 
on the assumption that all of the change in 7 during a merger 
can be attributed to the BHs, i.e., that cusp slopes remain 
unchanged the absence of BHs. This is known to be the 
case in equal-mass mergers between galaxies with p ower-law 
cusps (iBarnes 19991; iMilosavljevic & Merritt 200l|) , though 



in mergers with extreme mass ratios, features c an appear 
in the density profile that are not due to BHs (Merritt & 



Cruz 2001 ). We speculate that the break radii in some of the 



galaxies in our sample may be due to this process, particu- 
larly those galaxies (e.g., NGC 3640, 4168) where rb greatly 
exceeds the radius of gravitational influence of the BH. A'^- 
body simulations of cumulative mergers with unequal BH 
masses will be needed to assess this hypothesis. 

Our model presen ts an interesting contrast to that of 
van der Marel (1999| ), who proposed that cores (in the 
sense of constant-density regions) were present in all galax- 
ies ab initio, and that power-law cusps were generated by 
the growth of the BHs - roughly the opposite of our model 
in which BHs destroy pre-existing cusps. Van der Marel as- 
sumed that core mass correlated initially with bulge lumi- 
nosity as Mcoro ~ L^'^ and that M. oc L; hence Mcorc oc 
MI'^ , consistent with the correlation in Figure ^ if we iden- 
tify Mcoro with Mdef. We believe that this agreement is co- 
incidental. Van der Marel's model relates core mass to BH 
mass via an ad hoc postulate, while the model discussed here 
contains a mechanism for core formation. Van der Marel also 
ignored the effects of mergers. Nevertheless, van der Marel's 
model shows that our interpretation is not unique. 

If stellar cusps are destroyed by binary BHs, the same 
should be true of dark-matter cusps, like th ose predicted in 



CDM theories of structure formation (e.g., [Njaygxr o, Frenk 
fc White (19961 ); [Moore et al. (1998| ); pollock et al. (2001| )). 



Destruction of dark matter cusps could be very efficient if su- 
permassive BH s were present in d ark matter halos at larg e 
redshifts (e.g. [Fan et al. (200l|); [Haiman fc Loeb (2001^ ; 



Menou, Haiman & Narayanan (2001)) due to the cumula- 



tive effect mentioned above; furthermore, dark matter cusps 
would not be regenerated the way that stellar cusps might 
be via star formation. If our model for the formation of 
stellar cores is correct, we would predict tha cores of non- 
interacting CDM should be about as large as those observed 
in the stars, and perhaps larger. 

We thank L. Ferrarese for her generous help with the 
data reduction and interpretation. The image of M87 used 
for calculating the mass deficit of that galaxy (Section ^ 
was kindly provided by P. Cote and A. Jordan in advance 
of publication. This work was supported by NSF grant AST 
00-71099 and NASA grants NAG5-6037 and NAG5-9046 to 
DM. 
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15.0 


-18.9 


0.11 


1.69 


7.76 


8.22 


NGC 


4589 


22.0 


-20.0 


0.18 


1.05 


8.22 


8.78 


NGC 


5077 


34.0 


-20.4 


0.54 


1.15 


8.78 


9.51 


NGC 


5198 


34.1 


-20.0 


0.11 


0.92 


8.09 


8.63 


IN v_j V_/ 


OOUo 


26.6 


-19.7 


0.10 


1.84 


8.37 


8.01 


NGC 


5370 


41.3 


-19.0 


0.32 


1.56 


7.49 


8.63 


NGC 


5557 


42.5 


-21.2 


0.44 


0.87 


8.64 


9.62 


NGC 


5576 


25.5 


-20.3 


0.24 


1.37 


8.00 


9.19 


NGC 


5796 


36.5 


-20.7 


0.18 


1.18 


8.60 


9.15 


NGC 


5812 


26.9 


-20.3 


0.25 


1.71 


8.15 


8.83 


NGC 


5813 


32.2 


-21.1 


0.43 


0.31 


8.40 


9.66 


NGC 


5831 


27.2 


-19.9 


0.22 


1.42 


7.72 


8.67 


NGC 


5898 


29.1 


-20.4 


0.33 


1.34 


8.30 


9.20 


NGC 


5903 


33.9 


-20.9 


0.78 


0.77 


8.41 


9.57 


NGC 


5982 


39.3 


-20.9 


0.45 


0.49 


8.71 


9.58 


NGC 


6278 


37.1 


-19.7 


0.10 


1.53 


7.64 


8.59 


UGC 


4551 


23.6 


-18.7 


0.24 


1.36 


8.02 


8.87 



Table 1. Galaxies with 7inin < 2; D is distance in Mpc; rj, is 
break radius in kpc; 7inin is the minimum logarithmic slope; M, 
and M^^f are in solar masses. 



70 


E 


SO 


a = 4.8 


a = 3.75 


2.00 


0.93 ± 0.10 


-0.10 ±0.14 


0.91 ±0.09 


1.16 ±0.12 




(10.5 ± 1.8) 


(3.8 ± 0.86) 


(9.75 ± 1.78) 


(10.7 ±1.7) 


1.75 


1.07 ±0.15 


-0.78 ± 0.46 


0.86 ±0.24 


1.10 ±0.30 




(2.75 ± 0.78) 


(0.35 ±0.19) 


(2.17 ±0.50) 


(2.36 ±0.56) 


1.50 


1.56 ± 0.37 




1.58 ±0.35 


2.02 ± 0.45 




(0.49 ± 0.26) 




(0.46 ± 0.23) 


(0.55 ±0.25) 



Table 2. Linear regression fits to the (log A/, , log Afdcf ) data for 
three values of the fiducial logarithmic slope 70 . Values in paren- 
theses are M^^f/M, interpolated from the fit at M, = lO^AfQ. 
Fourth and fifth columns are, respectively, fits of the entire data 
set (inc luding the galaxies with dY namical estimates of A/^i^f) us - 



ing the Ferrarese & Merritt (2000) and the Gebhardt et al. (2000) 



values of the A/. — tr relation exponent a. 



